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I. INTRODUCTION 

The precise deterinination of system parameters by mea- 
surements is the basis of many applications in physics and 
beyond. Quantum mechanics offers a way to enhance mea- 
surement sensitivity by lowering the theoretically achievable 
limits. It is a long-standing goal in the quantum parameter 
estimation field to find feasible ways to reach this so-called 
Heisenberg limit. We call the number of physical systems 
used to measure the parameter (e.g., the number of atoms or 
photons) the resource size. In the Heisenberg limit, the mea- 
surement precision scales with the inverse resource size. In 
contrast, the classical shot-noise limit scales only with the in- 
verse square root of the resource size. Examples of emerg- 
ing quantum technology applications based on this enhance- 
ment are improved frequency standards [1], the construction 
of quantum clocks and their synchronization [2], the detection 
of weak forces [3], or the improved resolution and signal-to- 
noise ratio in image reconstruction [4, 5]. A straightforward 
proposal to achieve the Heisenberg limit in principle by using 
maximally entangled states faces the problem that these states 
are also more sensitive to background noise and hence deco- 
herence. In general this leads to a diminishing or complete 
cancellation of the overall improvement of these schemes [6]. 
In this paper, we present a scheme based on cluster states, 
which are a class of entangled states. They are well known 
in the context of quantum computation [7] but, to our knowl- 
edge, their relevance in parameter estimation theory has not 
yet been investigated. Cluster states prove very robust against 
many sources of noise [8, 9]. We combine this robustness 
with a scheme capable of estimating a phase parameter of the 
system evolution. 

The main results of our proposal can be summarized as (i) 
cluster states attain the Heisenberg limit for a setup with a 
three-body Hamiltonian and no decoherence. In atomic sys- 
tems, cluster states have been implemented with neutral atoms 
in an optical lattice [10]. Also many other proposals on how 
to create cluster states exist, for example, with ions [11], in 
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cavities [12, 13], in charge and flux qubits [14], and with the 
help of a linear quantum register [15, 16]. Proposals for im- 
plementing Hamiltonians of the form considered in this paper 
have been put forward for neutral atoms in optical lattices [17] 
and cold polar molecules [18]. It should be noted that re- 
cent works have shown that the Heisenberg limit is not the 
ultimate lower bound in quantum metrology, but that it can 
be beaten by implementing Hamiltonians with symmetric k- 
body interactions [19-21]. (ii) Afl required measurements are 
of a simple tensor product form involving only single Pauli 
measurements, (iii) Decoherence leads to a decrease in mea- 
surement precision but our scheme remains superior to various 
reference schemes. In implementations with cold atoms, the 
main contribution to decoherence is dephasing noise, which 
we study in detail in this paper. Furthermore, we consider 
depolarization and damping noise, which represent smaller 
contributions to the overall decoherence in typical realiza- 
tions. Our first reference system is given by a standard quan- 
tum metrology scheme with different Hamiltonian and initial 
states than the cluster state scheme. This will show that our 
scheme improves on the precision of these standard schemes. 
The second reference system consists of the Hamiltonian of 
the cluster state scheme but with different initial states. This 
will demonstrate that one can only expect this improvement 
for a suitable combination of Hamiltonian and initial state. 

Cluster states have been implemented in purely optical se- 
tups using polarization and momentum entanglement of pho- 
tons currently for up to six qubits [22-26]. The measurements 
required for our scheme can also be realized optically. We are 
not aware of proposals to implement the Hamiltonian of our 
proposed scheme in an optical system and developing such 
a proposal is beyond the scope of the present paper. Deco- 
herence in the optical case is mainly caused by photon loss, 
which is not discussed here [27]. 

This paper is organized as follows. After the introduction 
of parameter estimation theory and cluster states, we discuss 
in Sec. II how cluster states can be employed in parameter 
estimation schemes. It turns out that, in principle, it is pos- 
sible to achieve the Heisenberg limit with this cluster state 
scheme. We also define our reference systems. In Sec. Ill, 
we show how dephasing noise influences the precision of the 
measurements. We discuss an analytical model of the cluster 
state measurement scheme and present our main numerical re- 
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suits. These results are compared with an alternative estima- 
tion scheme and we find an overall improvement of precision 
even under noise. Two further noise channels, namely the de- 
polarizing and damping channels, are considered in Sec. IV, 
with qualitatively similar results. In Sec. V, we consider a 
further reference system subjected to all three noise channels. 
The scaling with the number of resources when we include 
noise is discussed in Sec. VI. There we show that our scheme 
offers advantages for small resource sizes. We conclude the 
paper in Sec. VII and explicitly derive an analytical solution 
for our model with resource size TV = 2 in an appendix. 



II. PARAMETER ESTIMATION SCHEMES 

In this section, we introduce the basics of parameter esti- 
mation theory and cluster states. Subsequently, we combine 
the two and introduce a scheme of parameter estimation with 
cluster states. Finally, we present the reference systems we 
will be using throughout the paper 



A. Parameter estimation theory 

We consider a system with Hamiltonian = x^o and 
density operator p{t) which evolves according to the von Neu- 
mann equation (h = 1) 

^it)^i[pit),H^]. (1) 

Parameter estimation aims to determine the value of the pa- 
rameter X by means of suitable measurements on many copies 
of the system. Repeated measurements of an operator O result 
in ly data points Oi. It is then possible to extract an estimate for 
X from this data by means of an estimator Xost(oi, . . . , o^). 
The uncertainty of the estimation is expressed as 

The denominator cancels superfluous "units" introduced by 
the estimator. If we assume an unbiased estimator, i.e., 
(Xost) = X' then the uncertainty, Eq. (2), reduces to the devia- 
tion of the estimator, which is a familiar measure of precision. 

If we let p{t) be a pure state, then it can be shown that the 
uncertainty 6x is bounded by [19, 28-30] 

'^^ - 

A Hq is the sta ndard deviation of Hq, i.e., AHq = 
\/ (Hq) — {Hq)'^. The bound Eq. (3) is also called quantum 
Cramer-Rao bound. There are two routes available to min- 
imizing the uncertainty. The factor in the denominator 
of Eq. (3) is the well-known statistical improvement with the 
number of independent runs i^. In general, this 1 / scaling 
is attained in the limit of many runs of the experiment [30]. 



Here we are not interested in this trivial improvement. The 
second and more interesting way is to maximize AiJo- 

Consider N replicas of a physical system evolving accord- 
ing to the Hamiltonian = xHq with 

N 

Ho = Y.hi'\ (4) 

j=i 

where h^^ is the Hamiltonian of the jth replica. We assume 
that each replica undergoes the same evolution so that each 
is of the same form. The maximal deviation of the com- 
posite Hamiltonian Hq is then attained by evolving the maxi- 
mally entangled state 




where \\min) and jAmax) are the eigenvectors of /iq with the 
smallest and largest eigenvalues Amin and Amax, respectively. 
This state yields AHq = A^(Amax — Amin)/2, and we see that 
the generalized uncertainty, Eq. (3), is given by 



max Aixiin ) 

Such a scaling with the inverse of the resource size, here N, 
is called Heisenberg scaling and the corresponding precision 
measurement is said to attain the Heisenberg limit. Note that 
this limit can be reached by using only separable measure- 
ments on each constituent [29]. In contrast, if we use uncor- 
rected states initially, the scaling is proportional to 1/ VN, 
which is known as the standard quantum limit. Quantum 
metrology offers an improvement of 1/\/N over classical 
schemes for Hamiltonians of the form in Eq. (4). However, 
in realistic setups the von Neumann equation is rarely exactly 
realized owing to the presence of environmental noise. Intu- 
itively one would expect noise to deteriorate the precision. In 
fact, Huelga et al. [6] and others [31, 32] found that decoher- 
ence rapidly destroys the improvement gained in a straightfor- 
ward implementation with maximally entangled states. They 
also point out the existence of partially entangled states which 
do improve the precision even under noise. This paper aims to 
improve on this result by introducing the more robust cluster 
states into quantum metrology. 

B. Cluster states 

Cluster states are a particular type of many-body entangled 
states [8, 33]. Their mathematical description is based on the 
notion of graphs. A graph G consists of a finite nonempty set 
V of N vertices together with a finite set E of m unordered 
pairs of distinct vertices from V [34]. Elements of E are also 
called edges since they join two vertices. Based on this notion 
we define the family of operators 
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where N{i) denotes the set of vertices sharing an edge with 

(i) 

the ith vertex, and a Pauli operator cri,j,,2 acts only on the ith 
vertex. A graph state |+g) is defined as the unique eigenstate 
with I+g) = I+g) for all operators K^^K The group Sg 
generated under multiplication by the set { K^"^") \ i £ V} is 
called the stabilizer of the graph state. It is instructive to iden- 
tify the vertices with a two-level physical system, which we 
will henceforth denote as qubit. With this notation, an equiva- 
lent way of describing graph states is given by first initiahsing 
every qubit in the superposition (l/\/2)(|0) + |1)) and then 
applying a controlled phase gate between each pair of qubits 
connected by an edge in the graph [33]. 

For simplicity, in this paper, we reduce the class of all valid 
graphs to linear graphs, which results in one-dimensional (ID) 
cluster states. Here, all vertices are connected to exactly two 
neighbors, which are pairwise distinct (except the first and last 
vertex, which are connected to one neighbor each). In phys- 
ical terms one could imagine a string of atoms, for example, 
each interacting with its next neighbors only. This leads to the 
family of stabilizers 



(5) 



for 1 < i < A^, where N is the number of vertices (or, equiv- 
alently, qubits). We use the convention that ai^^ = cri^^^' = 
1. Cluster states are known to be remarkably stable against 
a range of sources of noise [8, 9, 35, 36]. This stability can 
be traced to the fact that local decoherence on an individual 
qubit in the composite state only affects the small set of its 
neighbors. For cluster states, the number of neighbors does 
not depend on the total number of qubits but is a constant of 
the underlying graph. In this work we will exploit this feature 
to achieve a higher measurement precision even under the in- 
fluence of noise. 



C. Parameter estimation scheme with cluster states 

We have seen in Sec. II A that superpositions of states with 
maximally separated eigenvalues of Hq together with appro- 
priate measurements are sufficient to reach the Heisenberg 
limit. In Sec. II B we have introduced ID cluster states as 
eigenstates of a family of correlators A''*^ In analogy to the 
simple case of independent one-body Hamiltonians for a sys- 
tem of N qubits, we sum the correlators of all vertices and 
interpret this sum as a Hamiltonian 



H 



II 



N 

2 ^ 

i=i 



N 

2 / ^ Z ^ X ^ z 



(6) 



Here, x is the parameter we would Uke to estimate and j 
counts the qubits. In an atomic optical lattice implementa- 
tion of this Hamiltonian [18], the parameter x would sensi- 
tively depend on the atomic scattering properties. In a purely 
optical setup, measuring x could reveal the values of high- 
order terms in the susceptibility of the medium with high 
precision. Let I+at) be a ID cluster state of N qubits so 
that X^*) |+7v) = I+at) for all i. It is straightforward to 



-n) is an eigenstate 



see that the state |— tv) := Y\^=i '^^'^'^ 
with K^"^") \—n) = — \ —n) for all i. To prove this note that 
cri"''' commutes with K^'''^ for i ^ j but anticommutes for 
i ~ j [37]. This is analogous to the eigenstates |0) and |1) of 

the single particle Hamiltonians /iq*' = ctI*' in conventional 
quantum metrology schemes [6]. We use a superposition of 
eigenstates, cluster states, of the underlying Hamiltonian as a 
resource, namely. 



1 



^/2 



(I- 



(7) 



Here and in the following, an index c on a quantity denotes 
that this quantity belongs to the cluster state setup. Subse- 
quently, the system evolves under Hjj for a time t before it is 
measured. The time evolution leaves the system in the state 



1 / 

As measurement operator we employ 



N 

n 



r(j) 



The expectation value and deviation of AI^ are 

(M,) = cos{Nxt), 
AM, ^ \sm{Nxt)\. 



(8) 



(9) 
(10) 



We now consider v independent runs of the experiment. In or- 
der to calculate the uncertainty, Eq. (2), inherent in this setup, 
we define an estimator implicitly via 



COs(iVXestO: 



(11) 



i=i 



(1) 

where is the measurement operator for the jth indepen- 
dent run. This estimator is unbiased, which can be seen by 
expanding the right-hand side of Eq. (11) around the actual 
coupling X up to first order and taking the expectation value 
combined with Eq. (9). This reduces the uncertainty, Eq. (2), 
to the deviation 



SXo = AXest 



AMc 



1 



y/T^Nt\ sm{Nxt) I ^tN ■ 



Hence, our setup achieves the Heisenberg scaling 1 /N in the 
decoherence-free case. In this sense, the measurement op- 
erator Mf. we chose is optimal for this setup. The factor 
1 / y/V originates from statistical averaging over v measure- 
ments. The improvement over the shot-noise limit is 
the result of the TV-fold increase of the frequency in the ex- 
pectation value of Mc- 

As with all phase measurement schemes, the periodicity of 
the expectation value, Eq. (9), and estimator, Eq. (11), result in 
an ambiguity when globally determining the value of x with- 
out prior knowledge. However, we expect our setup to be a 
good candidate for local phase determination when the initial 
value of the parameter is known with a good accuracy [38]. 
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D. Reference systems 

In the following sections, we compare the attainable pre- 
cision of our cluster state setup with different reference sys- 
tems. The cluster state setup is defined by the Hamiltonian 
ff^ = Hi[, Eq. (6), the initial state |-0c), Eq. (7), and the 
measurement M^, Eq. (8). This setup is always kept invariant 
in these comparisons. 



III. PARAMETER ESTIMATION UNDER DEPHASING 
NOISE 

Physical systems are often exposed to various sources of 
noise. In this section, we focus on individual dephasing of 
the qubits. Dephasing is modeled by replacing the unitary 
evolution Eq. (1) with the master equation 



N 



(14) 



1. Reference system (1) 

The first reference system is the standard system used in 
parameter estimation. It is defined by the Hamiltonian 



N 



In this paper, we choose the Hamiltonians Hj or Hjj for 
depending on the reference system we study. The parame- 
ter 7 denotes the strength of the dephasing. To achieve the 
Cramer-Rao bound it is necessary to repeat the single mea- 
surements many times. As in [6], we assume a total run time 
T of the whole experiment with the duration t of each individ- 
ual realization. This results in the total number of experiments 
1/ = T/t. 



We then use two different initial states, which evolve under 
this Hamiltonian, and corresponding measurement operators. 
These are given by 



i=(|0)«^ + |l)«^ 



1 



|^„)=(^(|0) + |l)) 



N 

M^ = Y[<jy\ (12) 

N 

M, = "^ai'\ (13) 

3 = 1 



The indices m and u indicate quantities corresponding to max- 
imally entangled and uncorrelated states, respectively. The 
measurement operators Mm and Mu are optimal for attaining 
the lowest deviation in the decoherence-free case with the re- 
spective initial state [29]. If we include dephasing noise in this 
setup, we recover the result in [6] that maximally entangled 
states lead to the same minimal deviation 5x as uncorrelated 
states. The origin for this behavior is the A^-fold increase in 
the decoherence rate for maximally entangled state compared 
to uncorrelated states. 



2. Reference system (II) 

The Hamiltonian for our second reference system is given 
by Hjj, Eq. (6). As in reference system (I) we use maximally 
entangled and uncorrelated states as initial states for the pa- 
rameter estimation. This means that we can compare cluster 
state parameter estimation with parameter estimation with dif- 
ferent initial states under the same Hamiltonian Hjj. Also the 
measurement operators for the respective states are the same 
as before. Hence, reference system (II) is defined by Hamil- 
tonian Hjj, Eq. (6), together with Eqs. (12) and (13) as initial 
states and measurements. 



A. Analytical solution for N = 2 

In this section we discuss an analytical solution of the clus- 
ter state estimation problem with the master equation (14) for 
the case iV = 2. In the Appendix we present the derivation of 
this explicit solution. For higher TV our direct approach would 
also yield analytical solutions, however, the complexity of the 
calculation rapidly becomes intractable. For N ~ 2, the ex- 
pectation value and deviation of Mc are given by 



{Mc){t) =^ e-^* {cos{nt) + ^ siii(m) 



e^S/l 



-2-)t 



(cos(m) + ^sin(m))' 



^) siii(17t)| 



(15) 



(16) 



where S7 = \J Ax^ — 7^. The system frequency 57 depends on 
the dephasing rate but we can assume that 7/2% <C 1. This is 
justified by the fact that in parameter estimation, in general, it 
is desirable to have the parameter dominating over the noise. 
Comparing the deviation 5xc to the deviation of setup (I) with 
maximally entangled states. 



TtN\sin{Nxt)\ 



(17) 



with = 2, we see that the dephasing rate is reduced by a 
factor of 2. Furthermore, the additional term in {Mc){t) is of 
order 7/51, which is given by 



7 



2x 




Since 7/2% ^ 1, the additional terms in Eq. (16) are a small 
perturbation compared to the solution with maximally entan- 
gled states. Similarly, the frequencies in the reference sys- 
tem and the cluster state setup are the same to lowest order in 
7/2x- 
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We calculate the minimum of 5xc{t) by first finding the 
envelope function of Sxc{t)^ which is given by the solution of 
cos($) + asin($) = 0. We assume that a j/il and $ := 
ilt can be varied independently. Together with the equality 
cot(x — nnr) = cot(a;), where m is an integer, this yields 



o,.».-.(-»,.|+2,e,((2 



A: odd. (18) 



Plugging this into Eq. (16) results in the envelope of Sxc{t) 



^Xc,env(0 



2x/r(7 + AxH) 



(19) 



In the next step, we calculate the minimum of this envelope 
with respect to t, which is given by 




Since we have 7/2% <C 1, this minimal time to lowest order 
is given by tc,min ~ I/27 — 7/8x^1 which is approximately 2 
times the minimal time of setup (1) with maximally entangled 
states [6] if we further assume that I/27 S> 7/8x^- Hence, in 
the cluster state setup the evolution times can be longer than 
with maximally entangled states. It turns out that the minimal 
deviation (5xc,env(ic,min) is lower than the one for the standard 
scheme (1) with maximally entangled (or uncorrelated) states. 
If the condition, Eq. (18), is not satisfied exactly for a partic- 
ular set of parameters, then the actual lobes of the deviation 
still lie on the envelope but the minimum of these lobes does 
not coincide with the minimum of the envelope making the 
actual minimum of the deviation slightly larger. 

We are now able to calculate the relative improvement of 
this scheme with respect to the reference systems. The im- 
provement is quantified by 



e(7) = 1 - 



(5Xc,env(ic ,min I 



(20) 



ref 



where Sxrei is the minimum of the envelope of the reference 
deviation. With this definition an improvement of the cluster 
state scheme is given by positive values of e, whereas negative 
values indicate that the reference system performs better For 
both maximally entangled and uncorrelated states in setup (I) 
we have Sxief = <^Xm,min = SxuMn = ^/2jc/ NT. The ex- 
pansion of 5(7) in terms of the small parameter 7/x is given 

by 



e(7) 



1 




.X/ / 
(21) 

In the limit 7 = this series assumes an offset 1 - I/V2 w 
0.293. In contrast, in Sec. II we have seen that the cluster 
state setup and the conventional setup (I) yield the same result 
for 7 = 0. This discrepancy is owing to the fact that tc,min 
diverges as 7 ^ 0. For nonzero 7 we always find times for 
which Eq. (21) holds, although for small times the two devia- 
tions can lie arbitrarily close together. The position of the min- 
imum of the deviations increases with I/7, so as 7 decreases. 



the duration of the experiment quickly exceeds experimentally 
reasonable time scales. If we assume a fixed maximal evolu- 
tion time, the improvement can diverge from the constant for 
values of 7 smaller than 2 times the reciprocal of this evolu- 
tion time. 



B. Numerical results for N > 2 

In principle, the same analytical methods outlined above 
could be applied to > 2. Furthermore, the analytical re- 
sults could be compared to an analytical solution for reference 
setup (II), which would lead to qualitatively similar results. 
For simplicity, in the following we present numerical results 
of the cluster state setup and reference system (I) for higher N , 
and explain the findings with the insight we gained from the 
analytical treatment of the case N ~ 2. Our numerical results 
for N = 2 agree with the analytical solution above and known 
analytical solutions for arbitrary N with maximally entangled 
and uncorrelated states so we expect a good accuracy for the 
case N > 2. Similar to the case = 2 we find that the cluster 
state setup beats maximally entangled states in terms of low- 
est deviation even for > 2. In our calculations we limit the 
integration times to an A^-dependent time tf{N). This choice 
illustrates the effect of a finite evolution time in experiments. 

The results in Fig. 1(a) show the deviations for the A^ = 
7 cluster state setup compared to maximally entangled and 
uncorrelated states in reference setup (I). The envelopes of the 
deviation for maximally entangled and uncorrelated states are 
given by 



^Xu,o 



NVTt 



(22) 
(23) 



respectively. Note that both reference systems attain the same 
minimum at different times. The cluster state setup, however, 
attains a lower minimum than both, which is indicated by the 
black curve entering the shaded area in Fig. 1(a). The decay 
rate of the cluster state setup is approximately one-half the 
rate of maximally entangled states, which allows for a longer 
evolution time per experiment just as in the N ~ 2 case. 

The numerical analysis indicates that the envelope of the 
cluster state setup can be approximately expressed as 



I^Xc 



,Ar7t/2 



NVTt 



(24) 



This can also be extrapolated from the analytical envelope for 
N = 2, Eq. (19). To this end, in Eq. (19) we introduce the 
substitutions 2% Nx and 7 — ^ N^ /2. If we now assume 
^ /2Nx^ <^ 1, Eq. (24) follows. This approximation does not 
take into account the small oscillations we observe for 6xc 
around this envelope in Fig. 1(a). The minimum of Eq. (24) 
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FIG. 1: (Color online) Cluster states vs reference system (I). Plot (a) 
shows the rescaled deviation for the cluster state setup (solid line), 
and the envelopes for maximally entangled states, Eq. (22) (dashed 
line), uncorrected states, Eq. (23) (dashed-dotted line) in reference 
system (I), and cluster states, Eq. (24) (dotted line), all with N = 7. 
The shaded area is the area of possible improvement, which is pen- 
etrated by the cluster state scheme deviation. In (b) we plot the im- 
provement of the cluster state scheme vs maximally entangled states 
in reference system (I) for A'^ = 2, 3, 4, 5, 6, 7 (dark to light lines). 
Uncorrelated states lead to very similar curves. The dashed line in- 
dicates the analytical approximation, Eq. (21). The evolution time 
was kept constant for each at (TV) = [1/0.005A'^] so we expect 
major improvements for 7 > 0.005 in the shaded area. The horizon- 
tal dotted line marks the constant improvement 1 — l/\/2. For all 
curves x = 1 ™d in (a) 7 — 0.05. 



and the time at which the minimum is attained are given by 



A - / 

OXc.min — y rpj^' 

1 

^c.min — "TT • 

A'7 



(25) 
(26) 



This indicates that, for 7/x <C 1, we can expect the improve- 
ment, Eq. (20), of the minimal deviation to be e « 1 — 1 / V2 ~ 
29.3% compared to more common schemes with maximally 
entangled and uncorrelated states. This is consistent with our 
result for N ^ 2, where we obtained the same improvement 
to lowest order in j/x- Indeed if we plot the improvement for 
varying 7 for > 2, we recover this behavior for small 7/x 



(onset of the light shaded area in Fig. 1(b)). For Fig. 1(b) we 
fixed the evolution time to t — tf (N). This means that from a 
finite 7/ oc 1 /Nt f (TV) the minimum is expected to lie outside 
the integration domain, which leads to the observed deviation 
from the constant value in the white area (7 < 7/) of the fig- 
ure. For = 2 we know that the improvement will deviate 
from the constant value 1 — 1 /\/2 as 7/x ^ 1. 

We now focus on the substructure of kinks and humps of 
the = 2 curve (black line) in Fig. 1(b). Analyzing the 
position of the local maxima (humps) allows us to determine 
the frequency = V ~ T^- The humps are the result of 
the minimum of a lobe of the actual curve for 6x hitting the 
global minimum of the envelope of 5xc- This happens when 
the condition, Eq. (18), is satisfied at iniin — l/N-f. Hence for 
N ~ 2we rewrite Eq. (18) as 



2^ 



kn 

Y 



7 = 



2X 



(k odd). 



For the parameters in Fig. 1(b) the zeroth-order approxima- 
tion we used in this calculation locates the humps reasonably 
accurate for k > 3. For fc = 1 one would have to take into 
account higher-order terms in 7/O (cf. Eq. (18)). In between 
the humps for different k the kinks indicate where the lobe 
of Sx closest to the global minimum of the envelope is least 
optimal. We observe that qualitatively the same substructure 
persists for higher A^ but the humps slightly shift position and 
larger improvement sets in at lower 7. 



IV. DEPOLARIZATION AND PURE DAMPING 

In this section, we focus on two further noise channels, 
namely, depolarization and pure damping with the corre- 
sponding master equations 



dp 
dt 



N 



(t)=i[p,ffx] + f ^ 2a 



i=i 



'pa. 



(27) 



dt 



N 



3 = 1 



(28) 



respectively [35, 39]. Here we used cr± := (1/2) (cr^; ± ify)- 
In this section, we use system (I) as reference. Depolariza- 
tion occurs when the system interacts with a bath in the high- 
temperature limit T ^ 00. This drives the qubits into a com- 
pletely depolarized state, i.e., limt_,oo(o'z) — 0. Pure damp- 
ing, on the other hand, describes the decay of every qubit 
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FIG. 2: (Color online) Improvement of cluster states vs maximally 
entangled states in reference system (I) under depolarization. In this 
plot — 2, 4, 6 (dark to light lines) and x = 1- For uncorrelated 
states we get qualitatively similar results because both attain the same 
minimum. The shaded area indicates the region beyond the cutoff 
7/ = 0.005. 



into the state |0). Under the depolarizing channel, Eq. (27), 
both maximally entangled and uncorrelated states yield the 
same deviation as with dephasing (assuming the same mea- 
surements Mu and M^, respectively). This is the result of their 
measurements only extracting the "transversal" evolution de- 
pending on (jj,, whereas the depolarizing channel changes de- 
coherence in the longitudinal direction. On the other hand, 
the damping channel, Eq. (28), has the transversal effect of 
reducing the decay rate by one-half for maximally entangled 
and uncoiTclated states [35]. Hence for the polarizing channel 
we can use envelopes, Eqs. (22) and (23), and for the damping 
channel we replace 7 by 7/2 in these envelopes. This leads to 
the minima 
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TN 
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Both minima are the same and their value is, in fact, the same 
minimum we approximated for cluster states in the dephasing 
channel in the preceeding section. However, a priori it is un- 
clear how the cluster states will evolve under these different 
master equations. 

First we discuss numerical results for the depolarizing 
channel. In Fig. 2 we plot the improvement of cluster states 
vs maximally entangled states in reference system (I). Similar 
to dephasing noise the improvement saturates for intermedi- 
ate 7 before it starts to oscillate as the noise 7 approaches the 
strength of the parameter x- The constant improvement in the 
depolarization case is lower than for dephasing at e w 13%. 
We observe a very similar behavior if we take uncorrelated 
states as reference because their deviation attains the same 
minimum under depolarization as maximally entangled states. 
As expected, the improvement decreases for 7 below the cut- 
off 7/, which defines a maximal evolution time as before. 




FIG. 3: (Color online) Cluster state vs reference system (II) under 
dephasing noise, (a) shows the rescaled deviations and (b) the rela- 
tive improvement with uncorrelated states. The same parameters and 
notation as in Fig. I apply but here we expose maximally entangled 
and uncorrelated states to reference system (II). In contrast to Fig. 1, 
we do not plot envelopes in (a) but the actual deviation for the given 
parameters. 



For the damping channel we could not observe a global im- 
provement of the cluster state setup compared to either maxi- 
mally entangled or uncorrelated states in reference system (I). 
The deviation for cluster states is nearly identical to the one 
for maximally entangled states. This is caused by the reduced 
decoherence rate of both reference states under damping as 
mentioned above. Hence, in a system which is purely sub- 
ject to damping, cluster states offer no improvement in terms 
of minimal deviation over preparing a system in reference (I). 
However, experiments are often prone to a combination of er- 
rors and all other noise models we studied offered an overall 
improvement. 



V. RESULTS FOR REFERENCE SYTEM (II) 

In this section, we discuss all three noise models but now 
compare the cluster state setup with reference system (II). We 
start with the results for dephasing noise. In Fig. 3(a) we plot 
the deviations for maximally entangled states (dashed lines) 
and uncorrelated states (dashed-dotted) and compare it to Sxc 
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of the cluster states (solid). We do not use the envelope as 
before because we did not derive an analytical solution for the 
reference states in system (II). Instead the curves represent 
the actual numerical deviations for a specific set of parame- 
ters. Similarly to reference system (I), cluster states attain the 
lowest overall deviation. Maximally entangled states prove 
useless in this setup as their deviation is almost one order of 
magnitude larger than for cluster states or even uncorrected 
states. Their frequency also does not depend on TV, which 
is the fact that causes their advantage in the noise-free case 
of setup (I). The curve for uncorrected states exhibits a "bi- 
modal" structure, where the lobes follow two sets of envelopes 
with different dephasing rates. If we plot the improvement e 
in Fig. 3(b), it becomes apparent that the enhancement achiev- 
able with cluster states is even more pronounced than in ref- 
erence system (I). The nature of the sharp divide between the 
improvements for = 2, 3 and higher N is unclear. The rea- 
son for the larger improvement for 7 < 7/ = 0.005 in some 
curves is that the deviation of the uncorrected states reaches 
its minimum after the cluster states. In this regime of 7 the ref- 
erence system has not yet reached its minimum at the cutoff 
time tf, whereas the cluster state deviation is closer to its min- 
imum or has already reached it. Similar to the results with ref- 
erence system (I), we find plateaus for small 7, which indicate 
that we have again a constant offset as 7 approaches 0. For 
the resource sizes considered in the plot, these constant values 
correspond to improvements of approximately 40% and 60%. 

In contrast to the case with system (I), when we let the ref- 
erence states evolve under the depolarizing and pure damping 
master equations, we do observe a notable deviation improve- 
ment with cluster states under these noise models. The re- 
sults for damping in Figs. 4(a,b) suggest a steady increase of 
the improvement for higher particle numbers. One curiosity is 
the result for the N = 2 maximally entangled state in Fig. 4(b) 
(black curve). In this case, the expectation value of does 
not exhibit oscillations in time, which makes this state unsuit- 
able for the parameter estimation. For higher N we do recover 
oscillations but still see a large advantage of cluster states. The 
results for depolarization are similar with similar magnitudes 
of improvement for both reference states. 



VI. SCALING 





FIG. 4: (Color online) Improvement of cluster states vs (a) uncorre- 
lated states and (b) maximally entangled states in reference system 
(II) under damping. Both plots show A'^ — 2, 3, 4, 6 (from dark to 
light) with X = 1 and jf — 0.005 (marked by the shaded area). 



ation. We compare this expansion with the expansions of the 
deviation of maximally entangled and uncorrelated states in 
reference system (I), Eqs. (22) and (23), 




(30) 
(31) 



We have seen in Sec. II that, in principle, our scheme can 
achieve the Heisenberg limit in the number of qubits N for 
decoherence-free systems. For dephasing noise we have given 
an approximation of the envelope of the deviation Eq. (24), 
which results in a shot-noise-limited scaling, Eq. (25). In this 
section, we consider experiments with a fixed evolution time 
and derive the scaling properties of the deviation in this limit. 

We consider -jt to be a small parameter. Then we can ex- 
pand Eq. (24) to first order 



'^Xc.env — 



1 



NVTt 



0({N^t/2y 



(29) 



The first term on the right-hand side is a Heisenberg scaling 
with 1/iV and the second term is a constant offset of the devi- 



First we note that the constant offset (second term on the right- 
hand side) for cluster states in Eq. (29) is only one-half the 
offset for maximally entangled states in Eq. (30), and also 
smaller than the second term in Eq. (31). This guarantees that 
in this limit the deviation of the cluster states is always smaller 
than in both reference systems. Furthermore, the approxima- 
tion for cluster states is valid for {Njt/2)'^ <^ 1, whereas for 
maximally entangled states it is only valid for (N-ft)'^ ^ 1. 
For fixed -ft we can thus expect the deviation of the cluster 
states to scale with the Heisenberg limit (first term in the ex- 
pansion) for 2 times as many qubits as with maximally entan- 
gled states. In the limit considered here, uncorrelated states 
only scale with the shot-noise limit. For small numbers of 
qubits and small fixed jt the cluster state setup is superior in 
scaling to both reference schemes. 
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VII. CONCLUSION 

We have introduced cluster states into parameter estimation 
theory. By using "maximally entangled" cluster states as a 
probe in a system evolving under a three-body Hamiltonian, 
we could show that it is possible to achieve the Heisenberg 
limit of sensitivity to the parameter The addition of dephas- 
ing noise leads to decoherence in the system but our results 
indicate that the decoherence rate is reduced compared to the 
standard setup in quantum parameter estimation with maxi- 
mally entangled states. This leads to an improvement of the 
minimally achievable deviation. This improvement also per- 
sists when the cluster state setup is compared to different ref- 
erence setups. The results suggest that the improvement for a 
given number of qubits is almost constant when the strength 
of the noise is much smaller than the parameter, which is often 
desirable in experiments. Similar results hold when the qubits 
are subject to depolarization or damping noises. 

The creation of ID cluster states in optical lattices typically 
results in cluster sizes on the order of « 40 [10]. Also pro- 
posals for implementing three-body Hamiltonians have been 
devised for optical lattices [17, 18], which makes this sys- 
tem an ideal candidate for our measurement scheme. The 
achievable number of qubits could be increased even further 
by extending the present scheme to higher-dimensional clus- 
ter states. These can be created in optical lattices by using 
state-dependent lattice shifts. 

Analogous to parameter estimation with maximally entan- 
gled states, the Heisenberg scaling does not persist for large 
systems under decoherence. However, for small systems we 
have shown that the deviation of the cluster state setup scales 
at the Heisenberg limit. In the cluster state setup, this limit 
can be realized for 2 times as many qubits as with maximally 
entangled or uncorrected states. 
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ter Hamiltonian Eq. (6) for the case N = 2. Finding a so- 
lution for higher N is complicated by the fact that the basis 
of the underlying Hilbert space grows exponentially with N. 
In the following derivation we use the convention ctq •= 

where / is the identity. With this notation 



"■1.2,3 



^ x.y.z- 



the full Hamiltonian is given by 



X 



Hu^^i K 



'(1) 



X 

2 V "J 2 
The master equation takes the form 



.(1). 



.3 VW4'^+4^VW4- 



.(2), 



(2) 



2pit) 



(A.l) 



We expand p{t) e L(7Y®^) in the Pauli basis for two qubits 

3 

pit) = c^j{t)a^®aj, 

hj=0 

where (t) are the time-dependent coefficients of the corre- 
sponding basis vectors in L{TC^'^). For p{t) to be of unit trace 
we must have coo(i) = 1/2^ = 1/4. For the commutators 
in Eq. (A.l) we only need to consider terms in the expan- 
sion of p{t) where both Pauli matrices in the direct product 
do not commute with tri (g) (T3 and 0-3 ® ci . For the commu- 
tator with if'^-' (if'^') the commuting vectors are uo 'Si ctq, 
(To ® (T3, (Ti ® (To, and cti (g) (T3 (cto ® ctq, (To (g) (Ti , (T3 (g) (To, and 
(T3 (g) (Ti). Each commutator is then a Unear combination in the 
subspace orthogonal to the one spanned by the four commut- 
ing vectors. The incoherent part of the master equation can 
be easily evaluated by noting that cr'^^ p{t)a^^ flips the signs 

of the components which contain a[^^ or ctj'"'. These obser- 
vations allow us to rewrite the master equation as a system of 
coupled ODEs in the coefficients, which we write as a vector 
c{t) = [coi(t), co2{t), . . .]. Note that the component coo(i) is 
constant which does not affect the other coefficients so we do 
not include it in this calculation. In terms of this vector the 
system of ODEs becomes 



APPENDIX: ANALYTICAL SOLUTION OF THE 
DEPHASING MASTER EQUATION FOR = 2 



dc , , , , , 
-{t)=Ac{t). 



(A.2) 



To illustrate the behavior of the deviation Sxc, we derive an 
analytical solution of the master equation ( 14) with the clus- 



By evaluating the commutators and the incoherent part of the 
master equation we find for the matrix A, 
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This matrix is not time dependent so the solution of the 
ODE (A.2) is 

c{t) = e^'c{0), 

where c(0) are the components of the initial state. For the ID 
cluster state with N ^ 2 the initial values are given by 

COo(0)-i Cn(0) = -i C22(0) = ^, C33(0) = ^ 

and all other components vanish. We only need to deter- 
mine the quantity of interest {Mc)t = ti'[^4/o(0]' which can 
be derived from the solution c{t). The measurement opera- 
tor Mc = (73 (8i (73 leaves only one component of p{t) with 
nonzero trace 

(Afc)t = tr[c33(i)cro ® CTo] = 4C33(<) 



-7* 



(cos (nt) + ^ sin {nt) 



where = v4x 
the deviation 

SXcit) = 



2 — 7^. This leaves us with the solution for 

AA'/e 



/T/t\d{M,)t/dx\ 
2-^* -y/l - c-27t (cos(m) + ^ sm{nt)) ^ 



We note that for 7/0^0 this deviation approaches the de- 
viation for maximally entangled states with Hamiltonian Hj 
and a decoherence rate reduced by half [cf. Eq. (17)]. 



In a similar way one derives the solution for 
is given by 



1, which 



SXcit) 



-^^/y 1 - c--^* (cos(m/2) + ^ sin(m/2)) 
1+ jf^)sin(r!t/2)| 



^ cos{nt/2) - 



Comparing this form with Sxc (t) suggests an obvious exten- 
sion to iV > 2 by replacing the dephasing rate with Nj/2 
and the frequencies by Nil/2 in {Mc)t- However, we could 
not verify numerically that this generalized solution holds for 
N > 2. The numerical solutions for N > 2 follow this ex- 
trapolation closely but not exactly [cf. Fig. 1(a)]. 
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